Location of the Potts-critical end point in the frustrated Ising model 

on the square lattice 
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We report on Monte Carlo simulations for the two-dimensional frustrated J1-J2 Ising model on 
the square lattice. Recent analysis has shown that for the phase transition from the paramagnetic 
state to the antiferromagnetic coUinear state different phase-transition scenarios apply depending on 
the value of the frustration J2 / Ji . In particular a region with critical Ashkin- Teller- like behavior, 
i.e., a second-order phase transition with varying critical exponents, and a noncritical region with 
first-order indications were verified. However, the exact transition point [J2/Ji]c between both 
scenarios was under debate. In this paper we present Monte Carlo data which strengthens the 
conclusion of Jin et al. [Phys. Rev. Lett. 108, 045702 (2012)] that the transition point is at a value 
of J2/J1 ~ 0.67 and that double-peak structures in the energy histograms for larger values of J2/J1 
are unstable in a scaling analysis. 

PACS numbers: 64.60.De, 75.10.Hk, 05.70.Jk, 75.40.Mg 



I. INTRODUCTION 

The introduction of competing interactions in the clas- 
sical two-dimensional J1-J2 Ising model is accompanied 
by the appearance of new ground states and critical 
points at which the ground state shows a large degener- 
acy. Additionally ,the frustration can also affect the criti- 
cal behavior of the model. In particular, the emergence of 
varying critical exponents for the transition from a high- 
temperature paramagnetic phase to a collinear phase was 
observed numerically.^'^ 

However, a scenario of a continuous phase transition 
with varying exponents was in question after mean-field 
calculations by Moran-Lopez et al.^''^ gave evidence for 
a first-order transition in a certain regime of the frustra- 
tion. The first-order scenario was strengthened by Monte 
Carlo simulations which mainly focused on the evalua- 
tion of energy histograms in the vicinity of the transition 
temperature.^'^ 

Recently, we presented in Ref. 7 a full analysis for the 
phase transition using conformal field theory and exten- 
sive Monte Carlo simulations. The focus was on the 
derivation of the underlying field theory starting at the 
point of two decoupled Ising models at Ji =0, which 
is represented by a c = 1 field theory. By including 
the nearest-neighbor coupling Ji perturbatively, we ar- 
rived at an Ashkin- Teller field theory which has a central 
charge of c = 1 and is known to exhibit varying crit- 
ical exponents.® Moreover, the Potts-critical end point 
of this theory allows for the onset of a noncritical phase 
transition^' and could thereby explain the two different 
scenarios observed in the frustrated Ising model. The 
energy histograms showed two-peaked structures in the 
intermediate regime Ji/2 < J2 < Ji . This led to the con- 
clusion that the position of the critical end point was at 
Ji fa J2. Jin et al}^ agree with the general picture of two 
different transition scenarios and an underlying Ashkin- 



Teller field theory but argue that the first-order behavior 
is only valid up to J2 ^ 0.67 Ji. Using mainly argu- 
ments about the universality of the Binder cumulants, 
they show the equivalence of the point J2 ~ 0.67 Ji in the 
frustrated Ising model and the four-state Potts model, 
which marks the critical end point of the Ashkin- Teller 
field theory. ®"^° This universality for absolute values of 
the Binder cumulants was discussed extensively, e.g., for 
anisotropic couplings and different lattice systems in the 
Ising model and holds only under certain conditions. ^"^'^^ 



On the other hand, while the finite-size dependence of 
the histogram shape was analyzed for small parameters 
J2 = 0.6 Ji and 0.65 Ji,^'^ such an analysis was not per- 
formed for J2 > 0.7 Ji due to the increasing length scales 
needed to exhibit the double-peak structure in the first 
place. Thus, the onset of doubly peaked features in the 
histograms was interpreted as a signal for a first-order 
transition. In this work we present a finite-size analysis 
for the energy histograms on a larger scale and observe 
the vanishing of first-order signals for J2 > 0.67 Ji, in 
agreement with the conclusion of Jin et al.^^ that the 
first-order scenario ends at J2 ~ 0.67 Ji. Furthermore 
we analyze the trend of the critical exponents in the in- 
termediate regime 0.67 Ji < J2 < 1.2 Ji to verify the con- 
vergence of the Ashkin- Teller model and its exponents to 
the Potts-critical end point. 



After a short introduction to the model in Sec. II we 
present Monte Carlo results for histograms and critical 
exponents in Sec. Ill; in particular the scaling analysis 
will be presented as an example at the point J2 = 0.8 Ji. 
We summarize our results in the concluding Sec. IV. 
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FIG. 1. (Color online) Phase diagram for the frustrated Ising 
model for J2 > Ji/2. Two different regions are identified 
where the phase transitions show different behavior. In par- 
ticular the shaded area in the middle is discussed in this work 
by means of histograms and critical exponents {v, red trian- 
gles; 7, blue diamonds). Data for the critical temperatures 
(green circles) are from Ref. 7. 



FIG. 2. (Color online) Histograms for different lattice sizes 
recorded for critical (size dependent) temperatures at J2 = 
0.8 Ji. After the emergence of a double peak at intermedi- 
ate system sizes L = 1000, 1200 the vanishing of this first- 
order signature is observed for even larger systems (L = 
1500, 2000). 



II. MODEL 

The frustrated Ising model is described by the classical 
Hamiltonian 

-falsing = Jl ^ SiSj + J2 ^ SiSj (1) 
NN NNN 

which sums over all antiferromagnetic (Ji, J2 > 0) 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) 
interactions of spin variables iS^ = ±1 on a square lattice. 

Ground-state configurations arc given by a twofold de- 
generate Neel state for J2 < Ji/2 and a fourfold degen- 
erate coUinear state for J2 > Ji/2. At the critical point 
J2 = >/i/2 the transition temperature is suppressed to 
zero, and a ground-state manifold with a degeneracy of 
linear order is present.^ The phase transition to the Neel 
state belongs to the two-dimensional Ising universality 
class and will not be discussed any further. However, for 
J2 > two different phase-transition scenarios apply 
for small and large values of frustration Jnj J\'- Ashkin- 
TcUer like critical behavior with varying critical expo- 
nents for J2/J1 00 and first-order noncritical phase 
transitions for J2/J1 \ 0.5.^'^^ In Fig. 1 the critical tem- 
peratures are shown as green circles; the exact transition 
point between both regimes will be discussed in the fol- 
lowing section, and the area in question is marked in gray. 



III. MONTE CARLO RESULTS 

For the following numerical analysis we used a single- 
spin Metropolis update^^ with additional temperature 



exchange Monte Carlo steps^^ as it was also used in 
earlier works on the frustrated Ising model. 

To investigate the nature of the phase transition the 
first focus was on the computation of energy distribu- 
tions. In Ref. 7 we presented a histogram at J2 = 0.8 Ji 
for a lattice of linear size L = 1000 with a double-peak 
feature. This was interpreted as evidence for a first-order 
transition. For this work we computed systems with 
L = 500, 1000, 1200, 1500, 2000 and present in Fig. 2 
histograms at the size-dependent transition temperature. 
Surprisingly, the doubly peaked shape is observed only 
for intermediate system sizes 1000 < L < 1500. These 
first-order signatures appear at a particular crossover 
scale, but the distance between the two peaks decreases 
and vanishes completely when the linear system size is 
doubled again {L = 2000). A similar disappearance of 
a double-peak structure in energy histograms has been 
reported previously for the Baxter- Wu model. In that 
case, the two peaks approached each other in the ther- 
modynamic limit but were present for all finite system 
sizes. In contrast, the present situation where a double- 
peak structure emerges first at intermediate system sizes 
L « 1000 and then disappears again is highly unusual. 
Due to the vanishing double-peak structure the scenario 
of a first-order transition is no longer valid for the inter- 
mediate value of J2 = 0.8 Ji, and further analysis of the 
critical behavior is necessary. 

In light of the recent work by Jin et al}^ we checked 
the development of critical exponents for varying frus- 
tration J2/J1 \ 0.67 by analyzing our Monte Carlo data 
for system sizes from iV = 50 x 50 to = 500 x 500. 
A detailed investigation of correlation functions and the 
corresponding exponent 77 would not be meaningful since 
77 = 1/4 is expected to be constant in the Ashkin-Teller 
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FIG. 3. (Color online) Computation of critical exponents v 
and 7 from the Binder cumulants and susceptibility of the 
magnetic order parameter at J2 = 0.8 Ji. (a) and (c) show 
size-dependent values of the derivative [dU /dT]Tc and sus- 
ceptibility [x/Li'^]tc and algebraic fits that yield u = 0.78(1) 
and 7 = 1.36(1). In (b), additionally, a scaling collapse of the 
Binder cumulant is shown for given Tc and v; three original 
Binder cumulants are shown in the inset. 



model. By applying the scaling relation j/i/ = 2 — 1] the 
ratio of 7/1^ 7/4 is also fixed. However, the expo- 
nents 7 and v can be extracted separately. In particular 
the Binder cumulants^^'^-'^ Ub and the susceptibility x. 



for the collinear phase have been computed for this pur- 
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FIG. 4. (Color online) Flowgrams for the susceptibility x at 
the transition temperature Tc for several ratios Jij J\. Only 
for J2 > 0.7 Ji is the scaling nearly constant in L^"^^^ and 
Ashkin- Teller-like behavior obtained. 



pose. The order parameter is defined as 

m = m.j; + my , m^^y = j^^i-iy^'-^ Si 



(3) 



and satisfies (m) ~ for all temperatures in finite sys- 
tems. The results were double-checked for some param- 
eters with a slightly different definition of the order pa- 
rameter m2 = nii + mf, which was used in Ref. 11. In the 
scaling analysis for critical temperatures and exponents 
we did not notice any differences. 

The critical exponent ly of the correlation length ^ 
is calculated from the scaling of derivatives for Binder 
cumulants of different system sizes at the transition 
temperature:2'^ 



dU B / dT\Tc ~ aL 



(4) 



Hereby we neglect scaling corrections which have only a 
small infiuence on the resulting exponent, i.e., the de- 
viation is smaller than the fitting error. It should be 
noted that the analysis becomes more sensitive to cor- 
rections for J2/J1 — > 0.67. Nevertheless, the logarithmic 
corrections known to be present at the Potts point had no 
crucial influence in our analyses for J2/J1 > 0.7. As an 
example the scaling analysis is shown for J 2 = 0.8 Ji in 
Fig. 3(a): at the critical temperature Tc = 1.5679(4) Ji 
a critical exponent v = 0.78(1) is extracted. The same 
values are used in a finite-size-scaling collapse2'22^23 
the Binder cumulants vs L^l^ — Tc) in Fig. 3(b), and 
the good agreement in a wide temperature and lattice- 
size regime verifies the assumption of criticality in gen- 
eral and confirms the extracted values in particular. 
The same procedure was performed for several ratios of 
J1I Ji > 0.67, and all values of 1/ are shown in Fig. 1 
as red triangles. It is clearly observable that the expo- 
nents vary monotonically from the Ising value losing = 1 
for Ji = (right-hand side) towards the four-state Potts 
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value^^ I'Potts = 2/3 for J2/ Ji \ 0.67 (straight black ver- 
tical line) . A comparison with earlier works is ambiguous: 
Our values for v and J2 > 0.7 Ji arc in good agreement 
with the data presented in Ref. 24, although our values 
have smaller errors. However, more recent analyses at the 
point Ji = Ji yield a value of roughly v — 0.84(1)^^"^'' at 
Tc = 2.082(1) Ji, whereas our analysis yields a value of 
ly = 0.88(2) for a critical temperature Tc = 2.0819(4) Ji. 
This discrepancy shows the importance of the additional 
scaling-coUapse analysis which helps to refine the critical 
temperature and exponent. 

A similar scaling analysis is performed for the suscep- 
tibility and is also shown as an example at J2 = 0.8 Ji in 
Fig. 3(c). However, the critical exponent 7 also depends 
on the previously extracted v and is obtained without 
scaling corrections via 

xha^bL-y/'^. (5) 

The development of 7(J2/^i) is illustrated by the blue 
diamonds in Fig. 1. The limits are given by 7ising = 7/4, 
and 7Potts — 7/6.^^ Along the line of phase transitions 
the ratio of 7 and v is roughly constant at 7/4 which 
is valid for both the Ising model and the Ashkin- Teller 
model including the Potts-critical end point. This fact 
is also observable from the flowgram presented in Fig. 4: 
We present the flowgram of the susceptibility for differ- 
ent parameters J2/ Ji in addition to the scaling analysis. 
According to Refs. 22 and 28, the behavior of the flow of 
some observable connected to the phase transition should 
change significantly if the type of this transition is al- 
tered. Such an alteration is visible in Fig. 4 for the fiow 
at J2/J1 = 0.6, which is not constant in the chosen scal- 
ing over L~'^/^. For all other parameters, i.e., J2I J\ > 0.7 
the flow of the susceptibility shows only a small deviation 
from a constant scaling, which fits into the picture of the 
Ashkin- Teller critical behavior. 



IV. DISCUSSION 

We presented Monte Carlo data on the paramagnetic- 
coUinear phase transition in the frustrated Ising model. 
Our data are in agreement with the recent findings of 
Jin et al}^ and provide further support for the Ashkin- 
TcUcr nature of this phase transition in a specified pa- 
rameter region. In particular we analyzed the evolution 
of the critical exponents v and 7 from their Ising val- 
ues (at J2/J1 Z' 00) to the four-state-Potts values (at 
J2I J\ ~ 0.67). Thus, the critical frustration parameter 
of [J2/Ji]c ~ 0.67 proposed by Jin et al.^^ is consistent 
with our analysis. A comparison of the absolute values 
of v with estimates given by Landau and Binder^'* shows 
good agreement and smaller errors for our data. How- 
ever, at the present level of accuracy the values deviate 
from recent findings of several other groups. We also 
presented a finite-size analysis for the energy distribu- 
tion at J2 = 0.8 Ji and showed that first-order signals in 
the histograms vanish for linear system sizes L > 1200, 
and therefore the scenario of a critical phase transition is 
strengthened also by the absence of a latent heat. 
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